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Abstract. - The electric conductance of a strip of undoped graphene increases in the presence of 
a disorder potential, which is smooth on atomic scales. The phenomenon is attributed to impurity- 
assisted resonant tunneling of massless Dirac fermions. Employing the transfer matrix approach we 
demonstrate the resonant character of the conductivity enhancement in the presence of a single 
impurity. We also calculate the two-terminal conductivity for the model with one-dimensional 
fluctuations of disorder potential by a mapping onto a problem of Anderson localization. 



A monoatomic layer of graphite, or graphene, has been 
recently proven to exist in nature [1-4]. Low-energy ex- 
citations in graphene are described by the "relativistic" 
massless Dirac equation, which gives us theoretical insight 
into exotic transport properties observed in this material. 
Undoped graphene is a gapless semiconductor, or semi- 
metal, with vanishing density of states at the Fermi level. 
One of the first experiments [2] shows that the conduc- 
tivity of graphene at low temperatures takes on a nearly 
universal value of the order of 4e^//i and increases if a 
doping potential of any polarity is applied. Systematic 
dependence of the minimal conductivity on the sample 
size has been observed recently in Ref. [17]. 

The peculiar band-structure of the two-dimensional car- 
bon, which mainly explains many recent experimental ob- 
servations, has already been calculated in 1947 by Wal- 
lace [5]. Nevertheless, the universal value of the minimal 
conductivity is not entirely understood. Many recent the- 
oretical studies [8-13] address the problem of the finite 
conductivity of the undoped graphene by employing the 
Kubo formula. Other works [14-16] show that the conduc- 
tance G of a ballistic graphene sample (of the width W 
much larger than the length L) scales as G = aW/L with 
the coefficient a = Ae'^/wh. This value of a coincides with 
the prediction made for the dc conductivity of disordered 
graphene [6-9, 12, 13] and agrees with experiments done 
on small samples [17]. 

In this work we develop an extension of the transfer 
matrix formalism of Refs. [16, 18, 19] in order to include 
the effects of disorder. Our main result is the enhance- 
ment of the zero temperature conductance at low dop- 




Fig. 1: A ribbon of undoped graphene is contacted by two 
metallic leads. The charge carriers tunnel from one lead to 
another via multiple tunneling states formed in the graphene 
strip. For clean sample with L <^W the conductance G scales 
as G = goW/nL, where go = Ae^/h is the conductance quan- 
tum in graphene and L is the length of the strip. An impurity 
placed inside the strip enhances the conductance in a vicinity 
of the Dirac point e <C hv/L, provided the impurity strength 
is close to one of the resonant values. 



ing by an impurity potential, which is smooth on atomic 
scales. (Such potential corresponds to a diagonal term 
in the Dirac Hamiltonian [20]). Our results agree with re- 
cent numerical studies [21,22] and with a related work [23], 
where the conductivity enhancement by smooth disorder 
with infinite correlation range was predicted. 

We analytically calculate the two-terminal conductivity 
of a graphene sheet in a model with one-dimensional fiuc- 
tuations of the disorder potential taking advantage of a 
mapping onto a problem of Anderson localization. In this 
model the conductivity is found to increase as the square 
root of the system size. We find that quantum interference 
effects are responsible for the leading contribution to the 
conductivity of graphene near the Dirac point. 
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We start by considering the effects of a single impurity 
in the setup depicted in Fig. 1. At low doping the con- 
ductance is determined by quasiparticle tunneling, which 
is independent on the boundary conditions in y-direction 
if L <C W. (For illustrations wc choose periodic boundary 
conditions with W/L = 7). We find that a single impurity 
placed in an ideal sheet of undoped graphene modifies the 
tunneling states and leads to the conductivity enhance- 
ment provided the impurity strength is close to one of the 
multiple resonant values. Away from the Dirac point the 
presence of an impurity causes a suppression of the con- 
ductance. 

In this study we restrict ourselves to the single-valley 
Dirac equation for graphene, 



-ihv a ■ V* + = 



(1) 



where ^ is a spinor of wave amplitudes for two non- 
equivalent sites of the honeycomb lattice. The Fermi- 
energy e and the impurity potential V{x,y) in graphene 
sample (0 < a; < L) are considered to be much smaller 
than the Fcrmi-cncrgy Ep in the ideal metallic leads 
(x < and x > L). For zero doping the conductance 
is determined by the states at the Dirac point, e = 0. 
Transport properties at finite energies determine the con- 
ductance of doped graphene. 

The Dirac equation in the leads has a trivial solution 
^ ^ exp(±ifer) with the wave vector k = {kx,q) for the 
energy e = hv^k% + — Ep. In order to make our nota- 
tions more compact we let Uv = 1 in the rest of the paper. 
The units arc reinstated in the final results and in the fig- 
ures. For definiteness we choose periodic boundary condi- 
tions in y direction, hence the transversal momentum q is 
quantized as g„ — 2Trn/W, with n = 0, ±1, ±2, • • • ± M. 
The value of M is determined by the Fermi energy Ep in 
the leads, M = Int [W/A^], where Xp = 2it/Ep, and the 
number of propagating channels is given by A'' = 2M -|- 1. 

For £ -C Ep the conductance is dominated by modes 
with a small transversal momentum qn <^ kp. The corre- 
sponding scattering state for a quasiparticle injected from 
the left lead is given by 

m 

^i""^ = xY.^nme^'^-^ x>L, (2) 



where the summation extends from — M to M. The factor 
of 4 in the conductance quantum is due to the additional 
spin and valley degeneracies. 

In order to find tnm we have to solve the scattering 
problem. The solution becomes more transparent if one 
takes advantage of the unitary rotation in the isospin space 

= {(Tx + (Jz)/\/2, which transforms the spinors x and (f) 
into (1,0) and (0,1), correspondingly. We combine such 
a rotation with the Fourier transform in the transversal 
direction and arrange the spinors 



1 /-^ 



(5) 



in the vector '4>{x) of length 2N . Then, the evolution 

of tp{x) inside the graphene sample can be written as 
■j/>(x) = 7^i/>(0), where the 2A^ x 2A^ transfer matrix 7^ 
fulfills the flux conservation law Tja^T^ = Gz- In the 
chosen basis the transfer matrix of the whole sample is 
straightforwardly related to the matrices of transmission 
and reflection amplitudes. 



T = Tr 



(6) 



defined in the channel space. 

The equation for the transfer matrix follows from the 
Dirac equation (1), 



d% 
dx 



(7) 



where 5 is a diagonal matrix with entries qn and 1 is the 
unit matrix in the channel space. The elements of V{x) 
are given by 



^^m(*^) — 
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V{x,y). 



(8) 



For V{x) = 0, we denote the solution to Eq. (7) as 

Tj°) = exp [(ctx O g + iec^z O i) . (9) 

The matrix 7^"^ gives rise to the conductance of the bal- 
listic strip of graphene, which was calculated in Ref. [16], 
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(10) 



where fe„ = {y^k'j, - q^,q„) ~ {kp,qn), and 
1 / 1 



X = 



V2\ 1 



(3) 



The conductance of the graphene strip is expressed 
through the transmission amplitudes tnm in Eq. (2) by 
the Landauer formula. 
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50 = 4eV/i, 



(4) 



The Fermi-energy s in the graphene sample is a monotonic 
function of the doping potential. The zero-temperature 
conductance (10), plotted in Fig. 2 with the solid line, is 
minimal at £ = and corresponds to u = Ae^/irh. The 
minimal conductivity of ballistic graphene is due to the 
evanescent modes, which exponentially decay in the trans- 
port direction with rates qn- 

It is instructive to start with a simple impurity poten- 
tial, which is localized along a line x = xo, 



V{x,y) = a{y)5{x - Xq). 



(11) 
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Fig. 2: The conductance of the graphene strip, W/L = 7, with 
a potential interface (14) in the middle of the sample. The 
curves are calculated from Eqs. (4,6,12). The solid line shows 
the conductance of the ballistic sample, a = 0. The dashed 
line corresponds to the case of a resonant impurity strength, 
a = 7r/2 + nn. The averaged conductance for the stochastic 
model with a, which is uniformly distributed in the interval 
(0, 2-k), is plotted with the dotted line. 
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Fig. 3: The conductance of the graphene strip, W/L = 7, 
calculated from Eqs. (4,6,12) in the presence of a single delta- 
functional impurity (cf. Eqs. (11,16)) as a function of the im- 
purity strength. Different curves correspond to the different 
values of the chemical potential (Fermi energy) in the strip: 
eL/hv — for the lowest curve, eL/Tiv — 2 for the curve in the 
middle, eL/hv = 3 for the upper curve. The dotted lines are 
guides to the eye. 



In this case the transfer matrix of the sample reads 

r^Tl"X^e-^^'^^%^°l (12) 

At the Dirac point, e = 0, we find from Eq. (6) 

i^^ — coah{qxo)e''°' cosh{q{L ~ xq)) 

+ sinh{qxo)e''"sinh{q{L- xo)), (13) 

where the matrix elements of a are given by the Fourier 
transform (8). It is evident from Eq. (13) that the con- 
ductance at e = is not affected by any potential located 
at the edges of the sample Xq = or Xq — L. 

In order to maximize the effect of the impurity we 
let a;o = L/2 and calculate the conductance from 
Eqs. (4,6,12). We consider in detail two limiting cases for 
the y-dependence of a{y): a constant and a delta- function. 

For the constant potential, 

a{y) = a, (14) 
we have d = al, hence we find, at e = 0, 

G = 9oT.— J ■ 2 ■ (15) 

^ cos^ a cosh qnL + sm a 

Note that for any a the conductance at the Dirac point 
is equal, or exceeds its value for a = 0. Moreover, the 
conductance is enhanced to G = goN if the parameter a 
equals one of the special values a„ = 7r(n-|-l/2), where n is 
an integer number. This is a resonant enhancement, which 
takes place only in a close vicinity of e = 0. Taking the 
limit TV — > CX3 first, we obtain the logarithmic singularity 
at the Dirac point G = goiW/iTL) In |e|-|-0(l). The energy 



dependence of G for the resonant values, a = «„, is shown 
in Fig. 2 with the dashed line. 

For a stochastic model with a fluctuating parameter a, 
one finds a moderate enhancement of the averaged conduc- 
tance due to the contribution of resonant configurations 
with a ~ an- If the fluctuations of a have a large am- 
plitude (strong disorder), the contributions from different 
resonances are summed up leading to a universal result. In 
this case one can regard a as a random quantity, which is 
uniformly distributed in the interval (0, 2tt). The averaged 
conductance for this model is plotted in Fig. 2 with the 
dotted line. It acquires the minimal value G = gQ{W/2L) 
(cr = 2e^/ft.) at the Dirac point. Note that the averaged 
conductance is enhanced as compared to that of a ballistic 
sample for e < L/fiv. For large doping the situation is op- 
posite, i.e. the conductance is suppressed by the impurity 
potential. 

For the delta-function potential, 

a{v) = {^W/M)5{v-yo), (16) 

we find the elements of the matrix d in Eq. (12) as 

a„™ = (7/Af)e'(«"-'''")'^°. (17) 

Note that the ratio W/M = Xp remains constant in 
the limit — > oo. We calculate the conductance from 
Eqs. (4,6,12) and plot the results in Fig. 3 as a function 
of the parameter 7 for three different values of e. Since 
we assume periodic boundary conditions the conductance 
does not depend on the value of uq. 

We see that the effect of a single delta-functional im- 
purity is much smaller than that of the constant potential 
(14), but the main features remain. At zero doping, e = 0, 
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the conductance is enhanced for the special values of the 
impurity strength 7 = a„. Unlike in the previous case 
the height of the peaks is finite and is determined by the 
ratio W/L. (The effect is bigger if the impurity potential 
has a finite width.) It is clear from the lowest curve in 
Fig. 3 that in the stochastic model of fluctuating 7 the 
conductance at the Dirac point is necessarily enhanced. 

Away from the Dirac point the effect of the impurity 
is modified. The conductance G becomes an oscillating 
function of 7, and, for sL ^ 1, its value at finite 7 is 
smaller or eqiial the value at 7 = 0. 

We have to stress that the conductivity enhancement 
at the Dirac point induced by disorder potential relies 
upon the symplectic symmetry of the transfer matrix 
{—q)axTx{q) = Gx, which holds as far as the potential 
V{x, y) is a scalar in the isospin space. The corresponding 
microscopic potential is smooth on atomic scales even in 
the limit V{x,y) oc 5{i — rg). 

Let us illustrate the approach to disordered graphene, 
which follows from Eq. (7). The flux conservation justifies 
the parameterization 











cosh A sinh A 
sinh A cosh A 



U2 








. . (18) 
where [/'i,2, U{ 2 are some unitary matrices in the channel 

space and A is a diagonal matrix. The values A„ for x = L 
determine the conductance of the graphene strip 



1 



cosh X„ 



(19) 



The detailed analysis of Eq. (7) in the parameterization 
(18) is a complex task, which is beyond the scope of the 
present study. The problem is greatly simplified in the 
"one-dimensional" limit V{x,y) = V{x) due to the ab- 
sence of mode mixing. In this case the unitary matrices 
in the decomposition (18) are diagonal and the prime cor- 
responds to the complex conjugation. We parameterize 
Ui = diag {exp (z6'„)}, n = —M, . . . M, and reduce Eq. (7) 
to a pair of coupled equations for each mode 



dXn 
dx 

dx 



= g„cos26'„, (20) 
= £ - V{x) - q„ sin 20^ coth 2A„. (21) 



For V = s = 0, the transfer matrix fulfills an additional 
chiral symmetry, hence ^„ = 0. In this case the variables 
A„ grow with the maximal rate A„/.t = g„ as x increases, 
which corresponds to the minimal conductance. Any fi- 
nite doping, £ 7^ 0, or arbitrary potential V{x) violates 
the chiral symmetry and move the phases 0„ away from 
9n = 0. It follows from Eq. (20) that A„ < g„x, hence 
the conductance defined by Eq. (19) is enhanced above its 
value for ^ = e = 0. In general case of arbitrary V{x,y) 
the conductance is enhanced only on average, since a rare 



fluctuations with suppressed conductance become possi- 
ble. One illustration for the enhancement of the c;onduc- 
tance in the presence of mode mixing is provided by the 
lowest curve in Fig. 3. 

The; effcicts of individual impurities on the resistivity 
of graphene samples in strong magnetic fields have been 
demonstrated in recent experiments [24,25]. We, there- 
fore, believe that the phenomenon of the impurity-assisted 
tunneling considered above allows for an experimental 
test. 

Let us now give a brief analysis of th{^ c:onductivity in the 
model with a one- dimensional disorder, which is described 
by a white-noise correlator in the transport direction 

{V{x)V{x')) = ^^5{x-x'), {V{x)) = 0, (22) 

and is assumed to be constant in the transversal direction. 

Even though such a choice of disorder potential is clearly 
artificial, it gives rise to an an analytically tractable model. 
Due to the absence of mode mixing we can omit the in- 
dex n in Eqs. (20,21) and study the fluctuating variable 
A as a function of q and L. The two-terminal conductiv- 
ity a = LG/W, where G is found from Eq. (19), is given 
in the limit TV —> 00 by the integral over the transversal 
momentum 



goL 



f 

J — ( 



dq 



1 



27r cosh^ A(g,L) 



(23) 



We note that Eqs. (20,21) with the white noise po- 
tential (22) are analogous to the corresponding equations 
arising in the problem of Anderson localization on a one- 
dimensional lattice in a vicinity of the band center [26] . 

We look for the solution in the limit of large system 
size L ^ £ in which case the standard arguments can be 
applied. First of all, the variable A is self-averaging in 
the limit L ^ i, therefore the mean conductivity can be 
estimated by the substitution of the averaged value of A 
in Eq. (23), 



(A) 



Jo 



dx{cos29) ~ qL{cos2e), 



(24) 



where the mean value of cos 29 in the last expression is to 
be found from the stationary probability density P{9) of 
the phase variable. The main contribution to the integral 
in Eq. (23) comes from A ~ 1 since very small values 
of A are not affected by disorder. As the result we can 
let coth A ^ 1 in Eq. (21) and derive the Fokker-Planck 
equation on P{9) in the stationary limit i 



The solution to Eq. (25) has the form 



(25) 



POO 

P{e) QC dt exp [-As£t + Aql smt sin(i - 26^)] , (26) 
Jo 
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which leads to 



{\{q,L))=qL 



-4eH 



/i(4g^ sini) sint 



/o°°die-4^«/o(4g^sini) 



(27) 



where Iq, h stay for the Bessel functions. 

We notice that in the limit L ^ £ the integral in Eq. (23) 
is determined by the modes with qi <^1. For such modes 
we can let Ii{u) = u/2, Io{u) = 1 in Eq. (27) and obtain 



(A) = q'^lL 



1 



1 + (2£^)2 



(28) 



An interesting observation can be made at this stage. Ex- 
ploiting the analogy with Anderson localization a bit fur- 
ther we introduce a notion of the mode-dependent local- 
ization length ^ from the relation (A) = L/^. We, then, ar- 
rive at the standard result ^ = {e/q)'^£ (which means that 
the localization length is set up by the mean free path) 
only in the limit of large doping ei^ 1. On contrary, for 
e£ 1 we find the counterintuitive inverse dependence 
^ = {q^£)~^- This emphasizes once again an intimate re- 
lation of the underlying physics to the disorder-assisted 
timneling [27,28], which indeed suggests an enhancement 
of the length ^ with increasing disorder strength. 
Substitution of Eq. (28) to Eq. (23) yields 



o- = ego 




L~>£, 



(29) 



with the constant c ~ 0.303. Thus, the two- terminal con- 
ductivity in the model with one-dimensional fluctuations 
of the disorder potential increases with the system size 
without a saturation. The width of the conductivity min- 
imum is essentially broadened by disorder and is defined 
by the inverse mean free path %v/£ instead of the inverse 
system size fiv/L in the ballistic case. In the calculation 
presented above we have chosen to average A rather than 
cosh~^ A. This cannot affect the functional form of the 
result (29), however, the numerical constant c can slightly 
depend on the averaging procedure. 

Even though the localization effects are very important 
in the derivation of Eq. (29), the Anderson localization 
in its original sense is absent in the considered model. 
Indeed, Eq. (29) assumes that the conductance of the 
graphene strip decays as 

A generic random scalar potential in Eq. (1) would lead 
to the mode-mixing unlike the specific random potential 
(22) considered above. It is natural to expect, on the 
basis of the single-impurity analysis (16), that the mode- 
mixing will strongly suppress the effect of the conductivity 
enhancement. Nevertheless, the size dependent growth 
of the conductivity has been conjectured in Ref. [6] for 
a model of Dirac fermions in a generic two-dimensional 
scalar disordered potential. Moreover, the very recent nu- 
merical studies [29] provide a solid evidence of the loga- 
rithmic increase of the conductivity with the system size. 
The behaviour observed in [29] can be described by the 



expression of the type (29) provided \/L/£ is replaced by 
ln(L/£), where £ is inversely proportional to the potential 
strength. These results are in sharp contradiction with the 
recent work by Ostrovsky et al. [30], where a non- trivial 
renormalization-group flow with a novel fixed point corre- 
sponding to a scale-invariant conductivity a* ~ 0.6 is 
predicted. 

The conductivity enhancement discussed above re- 
lies upon the symplectic symmetry of the model (1). 
We should remind that quantum interference effects are 
also responsible for a size-dependence of conductivity 
of a disordered normal metal. In two dimensions, the 
zero temperature conductivity, which includes the weak- 
localization correction, takes the well-known form 



ne^Ts , 2e^ L 

0'2D metal = ± — T ' 

m nn 



(30) 



where n is the density of states at the Fermi level, Tg is 
the scattering time, and £ is the mean free path. The pos- 
itive sign in Eq. (30) corresponds to the case of a strong 
spin-orbit scattering [31] (the symplectic symmetry class). 
Thus, in a normal metal the symplectic symmetry also 
gives rise to the conductivity enhancement. This is in 
contrast to the orthogonal symmetry, which leads to the 
negative correction in Eq. (30). The derivation of Eq. (30) 
takes advantage of the small parameter {kF£)~^ and can- 
not be generalized to graphene at low doping. Neverthe- 
less the numerical results of Ref. [29] can be formally de- 
scribed by applying Eq. (30) beyond its validity range; i.e. 
in the situation when the Drude contribution (given by 
the first term in Eq. (30)) is disregarded in the vicinity 
of the Dirac point as compared to the weak localization 
term. We note, however, that the weak localization is not 
the only source of the logarithmic size dependence of con- 
ductivity in graphene [10]. 

In summary, an impurity potential, which is smooth 

on atomic scales, improves the conductance of undoped 
graphene. A confined potential can lead to a greater en- 
hancement of the conductance than the uniform doping 
potential. One single impurity can noticeably affect the 
conductance provided its strength is tuned to one of the 
multiple resonant values. We develop the transfer-matrix 
approach to the disordered graphene and calculate the 
two-terminal conductivity in the model of one-dimensional 
potential fluctuations. The resulting conductivity is fully 
determined by interference effects and increases as the 
square root of the system size. 
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through SFB 513. 



p-5 



M. Titov 



REFERENCES 

[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and 

A. A. Firsov, Science 306, 666 (2004). 
[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 

M. I. Katsnclson, I. V. Grigorieva, S. V. Dubonos, and 

A. A. Firsov, Nature 438, 197 (2005). 
[3] Y. Zhang, J. P. Small, M. E. S. Amori, and P. Kim, 

Phys. Rev. Lett. 94, 176803 (2005). 
[4] Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

438, 201 (2005). 
[5] P. R. Wallace, Phys. Rev. 71, 622 (1947). 
[6] A. W. W. Ludwig, M. P. A. Fisher, R. Shankar, G. Grin- 
stein, Phys. Rev. B 50, 7526 (1994). 
[7] E. V. Gorbar, V. P. Gusynin, V. A. Miransky, and 

I. A. Shovkovy, Phys. Rev. B 66, 045108 (2002). 
[8] N. M. R. Peres, F. Guinea, and A. H. Castro Neto, 

Phys. Rev. B 73, 125411 (2006). 
[9] K. Zicglcr, Phys. Rev. Lett. 97, 266802 (2006). 
[10] I. L. Aleiner and K. B. Efetov, Phys. Rev. Lett. 97, 236801 

(2006). 

[11] A. Altland, Phys. Rev. Lett. 97, 236802 (2006). 

[12] J. Cserti, Phys. Rev. B 74, 033405 (2007). 

[13] P. M. Ostrovsky, I. V. Gornyi, and A. D. MirUn, 

Phys. Rev. B 74, 235443 (2006). 
[14] S. Ryu, C. Mudry, A. Furusaki, A. W. W. Ludwig, 

arXiv:cond-mat/0610598 
[15] M. I. Katsnelson, Eur. Phys. J. B 51, 157-160 (2006). 
[16] J. Tworzydlo, B. Trauzettol, M. Titov, A. Rycorz, and 

C. W. J. Beenakker, Phys. Rev. Lett. 96, 246802 (2006). 
[17] F. Miao, S. Wijeratne, U. Coskun, Y. Zhang, C. N. Lau, 

arXiv:cond-mat/0703052. 
[18] M. Titov and C. W. J. Bccnakkcr, Phys. Rev. B 74, 

041401(R) (2006). 
[19] V. V. Cheianov and V. I. Fal'ko, Phys. Rev. B 74, 

041403(R) (2006). 
[20] N. H. Shon and T. Ando, J. Phys. Soc. Japan 67, 2421 

(1998). 

[21] A. Rycerz, J. Tworzydlo, and C. W. J. Beenakker, cond- 
mat/0612446. 

[22] J. A. Verges, F. Guinea, G. Chiappe, and E. Louis, 

Phys. Rev. B 75, 085440 (2007). 
[23] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98, 

076602 (2007). 

[24] E. H. Hwang, S. Adam, S. Das Sarma, A. K. Geim, 

arXiv:cond-mat /0610834. 
[25] F. Schcdin, K. S. Novoselov, S. V. Morozov, D. Jiang, 

E. H. Hill, P. Blake, A. K. Geim, arXivxond- 

mat/0610809. 

[26] H. Schomerus and M. Titov, Phys. Rev. B 67, 100201(R) 
(2003). 

[27] V. D. FreiUkher, B. A. Liansky, I. V. Yurkcvich, 

A. A. Maradudin, and A. R. McGurn, Phys. Rev. E 51, 

6301 (1995). 
[28] J. M. Luck, J. Phys. A 37, 259 (2004). 
[29] J. H. Bardarson, J. Tworzydlo, P. W. Brouwer, and 

C. W. J. Beenakker, arXiv:0705.0886. 
[30] P. M. Ostrovsky, I. V. Gornyi, and A. D. MirUn, 

arXiv:cond-mat/0702115. 
[31] S. Hikami, A. 1. Larkin, and Y. Nagaoka, 

Prog. Theor. Phys. 63, 707 (1980). 



p-6 



